Non-monotonic fluctuation spectra of membranes pinned 
or tethered discretely to a substrate 
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The thermal fluctuation spectrum of a fluid membrane coupled harmonically to a solid support 
by an array of tethers is calculated. For strong tethers, this spectrum exhibits non-monotonic, 
anisotropic behavior with a relative maximum at a wavelength about twice the tether distance. The 
root mean square displacement is evaluated to estimate typical membrane displacements. Possible 
applications cover pillar-supported or polymer-tethered membranes. 
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Introduction. — Thermal shape fluctuations of fluid 
membranes in the vicinity of a substrate depend both on 
the elasticity of the membrane and the specific type of 
interaction with the substrate For laterally homoge- 
neous substrates the combination of steric, van der Waals 
and electrostatic interactions determines the strength of 
these fluctuations. Such spectra, as measured experi- 
mentally using video microscopy 0, , decrease with in- 
creasing wave-vector since shorter wave-length fluctua- 
tions cost more bending energy. A qualitatively differ- 
ent type of interaction arises from tethering or pinning a 
membrane at discrete points to a substrate. By adjust- 
ing the length of the tethers the mean distance between 
substrate and membrane can be controlled. As tether- 
ing molecules mostly thiolipids are used which consist 
of a lipid taiLa hydrophilic spacer (e.g. peptides 0,0 
or polymers 0,0) an d a sulfur-based linker to the sub- 
strate. Likewise, end-functionalized membrane proteins 
can serve as tethers 

Micro- or nano-patterned substrates involving equi- 
distant silicon pillars or gold dots have so far mainly been 
used to study the interaction with cell membranes |9( or 
an actin cortex [Tol | but similar experiments with vesicles 
should become possible as well. If the membrane binds 
specifically and firmly to these structures, it is effectively 
pinned at these points. For future biotechnological ap- 
plications of such systems an understanding of how the 
thermal shape fluctuations are affected by tethering or 
pinning is of paramount interest. 

In this paper, we determine the spectrum of shape fluc- 
tuations of such membranes. Surprisingly, we find a large 
range of parameters for which the fluctuations exhibit a 
non-monotonic behaviour with a maximum at a wave- 
length of the order of the tethering or pinning distance. 
Previous theoretical work on membrane conformations 
on structured substrates focussed on groove-like or rough 
substrates 0, 0] or vesicles adhering strongly to chem- 
ically patterned substrates |l3|. In a formally related 
theoretical development, motivated by the plasma mem- 
brane of red blood cells, fluctuations of a compound sys- 
tem coupling a fluid membrane to the cytoskeleton were 



investigated 13 13 G3 

The model. - The membrane is linked at discrete 
points r a to a substrate by N springs as sketched in 
Fig. ^ The membrane surface is parameterized by a 
height profile h(r) = h(x, y) over a rectangular substrate 
of extensions L x and L y . The height h(r) refers to the 
displacement with respect to the rest length l of the 
tethers. The total energy functional 
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comprises membrane elasticity with bending rigidity k 
and harmonic tethers with strength K a . Assuming pe- 
riodic boundary conditions in the lateral direction, a 
Fourier expansion of the height profile reads 
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with k ee (k x , k y ) and k XiV = (0, ±1, ±2, . . .) • {2it/L x . v ) . 
Since h(r) is real, the complex Fourier coefficients 
obey h_ k 



- hi.. The spatial average of the height 
is real. The energy can be written as 



profile, h , 

E = | ^ k k , h w D kk ,h k , with a non-diagonal coupling 




FIG. 1: Sketch of a membrane tethered to a substrate at 
discrete sites. A membrane conformation is characterized by 
the height profile h(r), which denotes a displacement relative 
to the flat configuration at a height given by the rest length 
Zo of the tethers. 
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matrix _D k k , . Following a scheme introduced by Lin and 
Brown [l7| . a transformation to independent variables 
leads to 



E = c-Mc 



^^^^^ C^-i <p *p ' / 



(3) 



The components c,, (with r = 0,1,2,3,...) of 
the vector c can be grouped into three sectors, 
c= (g/iQ) {Re/iq}, {Im/i q }), where q runs through all 
independent, non- vanishing wave- vectors. With the def- 
initions 
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for r = , ,, 

4 t . n (4) 
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and 



s/2K c 



for r = (1 st sector) 
y/2K^ cos(q(r) • r a ) for the 2 nd sector (5) 
\]2K a sin(q(r) ■ r a ) for the 3 rd sector 

the matrix M can be simplified to take the form 
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can then be 
M _1 , where 



The fluctuation spectrum 
extracted from the inverse 

M r,r> = kh(°r C r') M 

Below we will compare the spectrum of the discretely 
tethered membrane with a simplified model called "con- 
tinuous springs" where a membrane fluctuates in a later- 
ally homogeneous harmonic potential of strength 7 with 
energy functional 



E^ = ~ Jdx [dy{n [V 2 h{v)\ 2 + 7 fr 2 (r)} 


The fluctuation spectrum of this system reads 
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The overall strength of the discrete and the continu- 
ous springs are comparable for 7 = L K ^ a _Q K a , 
since these parameters lead to the same effective spring 
constant of the spring brush as a whole. 

Fluctuation spectra. — We focus on a membrane at- 
tached to a quadratic array |2l| of equally strong (K a = 
K), equidistant tethers, with A as tether lattice constant, 
in the limit L xy — > 00. First, the case of pinning, i.e., 
infinitely strong springs {K — * 00) , is analyzed. Dimen- 
sional analysis yields the spectrum in the form 
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with the amplitude at the origin (k = 0) given by 
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FIG. 2: The scaling function g(k x A, k y A) of the fluctuation 
spectrum of a fluid membrane pinned by a quadratic array of 
infinitely strong springs. 



and a numerically determined prefactor s ~ 3.9 • 10~ 3 . 
The scaling function g(k x A, k y A) is shown in Fig.0 The 
fourfold symmetry of the spectrum reflects the symme- 
try of the tether array. This spectrum is non-monotonic 
with four relative maxima and four saddle points. The 
wavelength at the maxima is A max ~ 1.11 • 2 A . Since un- 
dulation modes with a wavelength twice the distance A 
fulfill a pinning condition (i. c. no displacement at the site 
of the tethers), a naive guess would yield A max = 2 A . For 
a quasi-one-dimensional system (with k y = 0) this would 
indeed be true |Tsj . However, for the full two-dimensional 
case, the interplay of all fluctuation modes apparent in 
the coupling matrix M © leads to deviations from this 
naive expectation towards a slightly larger wavelength. 

From this spectrum, one can derive various one- 
dimensional spectra shown in Fig. [3J Cuts through the 
origin and a maximum, and through the origin and a sad- 
dle point, respectively, show the non-monotonicity. The 
saddle point appears at a (by a factor of approximately 
\/2 ) larger wavelength than the maxima, corresponding 
to the larger tether distance in this direction. Further- 
more, the average of the 2-d spectrum with respect to the 
azimuth angle is displayed. In the averaged spectrum, 
the wavelength of the maximum is A;" ax ~ 1.2 • 2 A. It 
is larger than the wavelength corresponding to the rela- 
tive maximum in the 2-d spectrum but smaller than the 
wavelength corresponding to a saddle point. 

For finite spring constants K, the fluctuation spectrum 
of a membrane tethered by a quadratic array can be ex- 
pressed as 



l^ k l 2 ! 
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with a scaling function / that incorporates the spring 
elasticity. Since / is continuous and /(oo, k x A, k y A) = 
g(k x A, k y A), for sufficiently large spring constants K the 
non-monotonicity of the spectrum persists. In Fig.^J the 
wavelength-dependence of the azimuthally averaged scal- 
ing function / is shown for different values of K A 2 / k and 
compared to the spectrum (JHJl of a continuous confining 
harmonic potential. For weak springs, a tethered mem- 
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FIG. 3: One-dimensional plots of the universal 2-d fluctuation 
spectrum g(k x A, k y A) of a membrane pinned to a quadratic 
array of pinning sites (cf. Fig.|5Jl. Cuts in two main directions, 
through a maximum (curve 1) and a saddle point (curve 2), 
of the 2-d spectrum are shown and compared to an azimuthal 
average (curve 3). Curve 4 shows an azimuthal average of a 
spectrum for a hexagonal arrangement of pinning sites. The 
closer packed sites induce a six-fold symmetric 2-d spectrum 
situated below the spectrum of the quadratic system. The 
constant s for the hexagonal pattern, s( hcx ) ~ 2.5T0 -3 , is dif- 
ferent than in the quadratic case Jim . For comparability, the 
hexagonal spectrum is divided by the prefactor ( \ho\ 2 } p i nnc d 
of the quadratic case JSJ. 



brane behaves like a membrane with continuous confine- 
ment. For moderately strong spring constants, deviations 
from the continuous spectrum occur, with the spectrum 
of the discrete tethering being systematically larger than 
the continuous one. The positions of the relative max- 
ima in the 2-d spectra depend on KA 2 /k. For infinitely 
strong springs the maximum of the spectrum reaches its 
lowest possible wavelength. For decreasing spring con- 
stants the wavelength corresponding to the relative max- 
ima increases. Finally, below a critical value of the spring 
constant, KA 2 / 'k\ — 100, the relative maxima disap- 
pear and the spectrum decays monotonically. 

Real-space fluctuations. — We now determine the av- 
erage width of the membrane in real space. The height 
profile can be written as 



h{r) = 2w(r) • c 



(12) 



with [13 



w(r) = (l,{cos(q-r)},{-sin(q T)}) . (13) 

Given the inverse matrix M _1 , the mean square displace- 
ment reads 

(/i 2 (r)) = 4^iy r (r)<c r c r ,)«v(r) 

= 2k B T ^ti) r (r)M^ w r ,(r) . (14) 

r,r f 

Typical membrane elongations are represented by the 
root mean square displacement (RMSD) -J (h 2 (r)). Only 
if the maximum of this RMSD is significantly smaller 
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(1) discrete springs, A'A 2 /k— -oo 

(2) conti. springs, aa 2 /k = iooo 

(3) discrete springs, aa 2 /k = 100 

(4) conti. springs, a^A 2 /k = ioo 

(5) discrete & conti., aa 2 /* = io 

(6) free membrane, ka 2 /k = q 
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FIG. 4: The scaling function f(K A 2 /k, k x A, k y A) for the 
fluctuation spectrum of a quadratically tethered membrane 
for different spring stiffnesses 1111 . Azimuthal averages of 
the 2-d spectra are compared to results from the continuous 
spring model (|HJ. For comparision, the spectra for the con- 
tinuous model are scaled by the prefactor ( |/io| 2 ) pinned °f the 
discrete case. Curve 1 shows the azimuthal average of the 
2-d spectrum for pinning, i. e. infinitely strong springs. In the 
comparable continuous case (7 — > 00), the spectrum vanishes. 
Curve 2 depicts the spectrum for strong continuous springs: 
While the "discrete spectrum" is still undistinguishable from 
curve 1, the "continuous spectrum" is located far underneath. 
Curves 3 and 4 give the corresponding discrete and continu- 
ous spectrum for moderately strong springs, respectively. For 
weaker springs the two cases yield nearly the same data shown 
in curve 5. The spectrum of a free membrane without tethers 
is shown as curve 6. 



than the rest length of the springs, this model can be 
trusted quantitatively since it does not yet include the 
steric hindrance by the substrate. The maximum of the 
RMSD is a measure to estimate how far above a sub- 
strate the membrane has to be placed in order to avoid 
undesirable contact or adhesion of the bilayer to the sub- 
strate. 

For pinning (K — > 00), we find 



pinned 



kBT (x y\ 

- a ' u Ia'aJ (15) 



with a scaling function u(x/ A, y/ A) shown in Fig.|3J The 
maximum of the RMSD in the pinned case occurs in the 
center of the array and is given by 



pinned 
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with a constant p ~ 0.124 . 

For finite spring constants K, the maximum of the 
RMSD takes the form 



max 



{VWW} = p 



k B T ( KA' 
A • w 



V K 



(17) 



with a scaling function w. For weak springs, w(K A 2 / k) 
approaches ( ^j ) 1 ^, which follows from the 
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FIG. 5: Scaling function u(x/A,y/A) for the RMSD, 
v/ ( h 2 (r) ) , of a membrane pinned at a quadratic array of 
sites. 

known result 

(h 2 (r)) M = -M= (18) 

for a membrane bound in a harmonic potential |19| . 
For strong springs, w approaches 1. The cross-over be- 
tween these two scaling limits occurs around KA 2 /k — 
(pV8)~ 4 — 66.1, where the two asymptotes intersect. 

Summarizing perspective. — Discrete tethering of a 
membrane to a substrate has a profound implication on 
the fluctuation spectrum. For strong enough tethers, this 
spectrum becomes non-monotonic with a maximum de- 



termined by the spacing of the tethers. An interpretation 
of such a spectrum in terms of a continuous model would 
require a term — |<t|(V/i) 2 in implying a negative "sur- 
face tension" . For sufficiently weak tethers, however, the 
discrete tethering can indeed be replaced by a continuous 
harmonic confining potential. The explicit introduction 
of an additional regular (positive) surface tension, e. g. 
caused by the area constraint of a vesicle whose bound 
part is considered j2jj, into the energy JQl of our model 
poses no problems. 

Our quantitative data on the mean displacements 
caused by these fluctuations should provide valuable 
hints for the experimentalists, how large the rest length l Q 
of the pillars and tethers has to be in order to avoid con- 
tact with the substrate. For a smaller Iq, one should in- 
clude a direct interaction with the substrate which could 
be attractive due to van-der-Waals interaction and/or 
repulsive due to steric interactions. In either case one 
could include an effective potential V(h) to the energy 
functional. Its minimization will then lead to a laterally 
inhomogeneous "ground state" profile ^ ( r ) as m studies 
of membrane adhesion to structured substrates [l^ . In 
another extension of our model one can study the fluc- 
tuation of two almost parallel membranes connected by 
polymeric linkers. The fluctuation of the relative dis- 
tance between the membranes, i. e., the peristaltic mode, 
is then governed by an effective Hamiltonian Q and will 
show the non-monotonic behavior for strong enough link- 
ers as well. 
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